Spreading of a granular droplet 



Ivan Sanchez, 1 Franck Raynaud, 2 ^ Jose Lanuza, 2 Bruno Andreotti, 2 Eric Clement, 2 and Igor S. Aranson 3 

1 Departamento de Fisica, Universidad Simon Bolivar, Apartado 89000, Caracas 1080- A, Venezuela 
2 PMMH, UMR7636 (CNRS), ESPCI Univ. P6-P7, 10 Rue Vauquelm, 75005 Pans, France. 
3 Materials Science Division, Argonne National Laboratory, Argonne, IL60439, USA 

(Dated: February 3, 2008) 

The influence of controlled vibrations on the granular rheology is investigated in a specifically 
designed experiment. We study experimentally a thin granular film spreading under the action of 
horizontal vibrations. A nonlinear diffusion equation is derived theoretically that describes the evo- 
lution of the deposit shape. A self-similar parabolic shape (the "granular droplet") and a spreading 
dynamics are predicted that both agree quantitatively with the experimental results. The theoretical 
analysis is used to extract effective dynamic friction coefficients between the base and the granu- 
lar layer under sustained and controlled vibrations. We derive an empirical friction law involving 
external driving parameters and the layer height to represent our data. 

PACS numbers: 47.10.+g, 68.08.-p, 68.08.Bc 



Granular matter is a model material presenting the 
phenomenology of a wide class of complex fluids with 
a yield stress governing the transition between solid-like 
and liquid-like behaviors. In this respect, the most im- 
portant challenge is to establish microscopically founded 
constitutive relations describing the different phases as 
well as the features of phase transition itself. In the 
last decade, various important propositions were made 
to draw analogies with glassy states of matter or out of 
equilibrium thermodynamics [TJ |2] but a full understand- 
ing still remains an open question. An important reason 
has been the failure of classical rheometers, which do 
not lead to homogeneous flows [3]. For dense granular 
flows, the situation has progressed by the introduction of 
new geometries ([4] and references therein) like the in- 
clined plane configuration, in which the stress tensor is 
controlled. Now, it is known that sufficiently far from 
the jamming transition, the dense granular flow rheol- 
ogy is local in first approximation: it is governed by an 
effective friction coefficient increasing with the shearing 
rate, properly rescaled by the local confining pressure. 
The kinetic theory, valid in the gaseous regime, fails to 
describe this situation characterized by the existence of 
long term contacts and force networks. In particular, bi- 
nary collisions lead to a decrease of friction jS] instead 
of the observed increase. The most important open is- 
sue concerns the jamming/unjamming transition and the 
lack of identification of parameters controlling this hys- 
teretic transition [B] . It has recently been suggested that 
the dynamics in the metastable region could be dom- 
inated by elementary rearrangements, interacting non- 
locally [7] by fluctuations or by the coupling with the 
modes of vibration [8^ (from internal or external origin). 

Here we have designed a specific experiment to mea- 
sure the influence of controlled vibrations on the granular 
rheology. We hereafter study the unjamming induced by 
horizontal shaking due to the transfer of energy by fric- 
tion. Note that the majority of previous experiments 




FIG. 1: (a) Sketch of the experimental setup, (b) Laser profile 
as a function of time: the top one corresponds to the granular 
film prior to shaking, the three below are subsequent deposit 
shapes after the vibration onset. Due to the laser inclination, 
the height is magnified by a factor ~ 10. (c) Forces acting on 
a vertical slice of the drop, (d) Analogy with a solid slider on 
a vibrating tray. 

have used vertical vibration to unjam and agitate granu- 
lar matter but relatively few have used horizontal shear 
H EH El E2], see also [T3] for review. We focus on 
two different questions. How does a granular film spread 
under the action of horizontal vibrations? How can the 
rheology be deduced from the spreading dynamics? 

Setup - The experimental apparatus is an horizontally- 
shaken tray (Fig.[T^) whose motion is a sinusoid of ampli- 
tude a and angular frequency to = 2nf. Driving frequen- 
cies / were between 15 Hz and 30 Hz and the maximum 
amplitude that we could reach is 1 mm. The substrate 
is a sand-cast aluminum plate with roughness of about 
10 (im. The tray is stabilized by 4 ball-bearing rings glid- 
ing on two rails. The tray is leveled horizontally within 
1/100°. The CCD video camera is moving along with 
the tray. The granular material is d ~ 300 /im diameter 
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faceted Fontainebleau sand of static and dynamical fric- 
tion coefficients p s = 0.66 and pd = 0.60. The following 
protocol is used to prepare a layer with reproducible ini- 
tial conditions: an initial mass of sand is poured in a rect- 
angular bottomless confining box fixed to the substrate. 
The box inner base defines the initial width L x — 40 mm 
in the vibration direction x and L y = 150 mm in the 
lateral dimension. To level off the deposit, the granular 
material is shaken vigorously for 1 min at 40 Hz. After 
removal of the confining box, an initial roof-like shape 
sand layer is obtained (Fig. [lj>, top). To keep a lateral 
confinement of the deposit, two small guiding sidewalls 
are fixed on the tray. The shape of the deposit is mon- 
itored by a laser sheet shined at a small angle in the 
central part of the layer. On Fig. [IJ>, we display below 
the initial profile, three subsequent laser slices observed 
for "vigorous" shaking. One sees that the granular layer 
looses its stability and consequently, spreads horizontally 
along the x direction. 

An inertial tribometer - From direct visualization of 
the layer motion using a fast camera (500 Hz), we have 
noticed that within an oscillation period, all the surface 
grains bear almost the same relative phase with respect to 
the plate. This suggests a close analogy with the motion 
of a solid slider. If a constant Coulomb static friction 
p s is assumed, the onset of relative motion between the 
substrate and the solid would correspond to a rescaled 
acceleration T = auj 2 / g equals p s . Actually, the situation 
turns out to be more complex. We leave a systematic 
study of the motion onset to a future report and we focus 
here on the high T regime for which the film is set into 
motion and eventually leaves the limit of the tray before 
stoppage. 

Let us consider the dynamics of a vertical slice of the 
granular droplet of width dx, of length L y and of height 
h{x,t) (Fig. [lp). Along the vertical direction, pressure 
P balances gravity g and reads: P = pg(h — z), where 
p is the density of the material. Assuming normal stress 
isotropy, the pressure gradient induces a driving force 
—L y dxpghd x h on the slice. The tray exert on it a re- 
sistive force Lydxr (r is the shear component of stress) 
that opposes, on the average, the spreading motion. To 
assess the momentum transfer due to the complex sliding 
dynamics between the droplet and the base, we introduce 
a friction coefficient p relating r to P. It is worth noting 
that p can depend on the basal pressure P, on h, on a and 
ui. The strong hypothesis is that it weakly depends on 
the relative velocity between the grains and the tray. We 
checked for instance that the introduction of a hysteresis 
between static and dynamic regimes does not affect sig- 
nificantly the results in the range of T investigated here. 

Then, the problem is completely similar to that of a 
solid block on an oscillating tray, submitted to a driving 
force F and a normal force N = mg (Fig. [TJi) . The anal- 
ogy is established through the dimensionless parameter 
F/N = —d x h. Our setup is in some sense the equivalent 
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FIG. 2: Spreading dynamics at a frequency / = 28 Hz for 
a section S = 90 ± 4 mm 2 , (a) Rescaled droplet shapes at 
different times t for vibration amplitudes a — 0.25 mm and 
a = 0.39 mm. (b) Fit residuals, (c) W(tf- W{0f as function 
of time t for various vibration amplitudes a. 



of the 'inertial tribometer' [2] designed in the context of 
solid on solid friction studies. We make the equations of 
motion dimensionless using a;" 1 as a characteristic time 
and a as a characteristic length. In the sliding regime, 
the equation governing the evolution of the dimensionless 
velocity difference v between the slider and the tray is: 

v F 



Tv = sin t 



/'T 



N 



(1) 



Blockage of motion occurs for v = 0. From rest (in the 
moving reference frame), motion starts when: |sint| = 
p/T. We have obtained exact solutions of the problem at 
the linear order in F/N. For the purpose of the present 
Letter, we only need the average sliding velocity v, which 
is written under the form: 

, = mal (2) 

p N w 

Simple calculations show that depending on the value of 
r//i, the block motion may undergo two different dynam- 
ical regimes: a continuously sliding at large acceleration, 
for which 
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and stop/start motions at smaller acceleration, for which 
p((f) — arcsin(/x/r)) 2 



P = 



2ttT 



(4) 
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where <f> is the root of the transcendental equation: 
cos(» - y/l - {n/T) 2 + n/T{<j> - arcsin(^/r)) = 0. As 
shown on Fig. |3ja) , (3 is an increasing function of T/fi 
that reflects the fraction of time during which the solid 
block slides. 

Spreading dynamics - A small aspect ratio H/W 
droplet allows for a depth integrated (Saint- Venant) de- 
scription. Assuming no significant variation of density 
p, the flux of grains across a vertical section of the drop 
can be approximated by auvh. Using Eq. (J2), the mass 
conservation law yields a nonlinear diffusion equation: 



d t h = d x (Ud x h 2 ) with U 



au>0(T / (i) 
2m 



(5) 



Note that the diffusion coefficient U is velocity indepen- 
dent. Assuming that U does not depend on height h, 
Eq. (|5| admits an exact parabolic self-similar solution: 



h(x, t) 
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with a spreading dynamics that can be written: 

W{tf = W{0f + 72SUt (7) 

It can be shown that sufficiently localized initial condi- 
tions converges at long time toward the self-similar so- 
lution - just like the convergence toward gaussians in 
linear diffusion. Experimentally, the profile shape is ex- 
tracted from images of the laser trace, using a standard 
correlation technique. We have fitted it by a second or- 
der polynomial in order to determine the width W(t) 
and the maximal height H(t). Fig. [2] shows the pro- 
files measured at different times, rescaled by W and H. 
Remarkably, all experimental points collapse on the pre- 
dicted parabolic shape (Eq. [6]) . Residuals of the differ- 
ence between the fitted curve and the experimental data 
shows a deviation to a parabola of less than 4% in the 
central part and around 10% on the edges, for all data 
presented here. As expected from mass conservation, the 
cross section area remains constant in time and is equal 
toS = 2H(t)W{t)/3. 

Several curves W(t) 3 — IU(0) 3 obtained for / = 26 Hz 
are displayed on Fig.|2p. They present a quasi linear time 
dependence whose slope increases with the vibration am- 
plitude a. We have chosen here initial width values W(0) 
such as to avoid the non-universal initial transient. Note 
that in our experiment, because of the lateral boundaries, 
we have evidenced weak but systematic transverse cur- 
vatures: the part in the center being the fastest spread- 
ing one. Strictly speaking, the spreading dynamics is 
three-dimensional but may be approximated by a two- 
dimensional description since the curvature effect is weak. 
So, theoretical predictions of the shape (Eq. [6]) and the 
spreading law (Eq. [7]) are both in excellent agreement 
with experimental findings, Fig. [2] 



Parametric study - Using Eqs. (|5| and the slope 
of the relation between W 3 and time t can be written as: 



dW 3 _ 18gST 
M 



dt 



(8) 



From this relation we get an effective value of the friction 
coefficient \x. Therefore, as initially desired, this spread- 
ing experiment under controlled vibrations can be used 
as a "granular tribometer" in order to monitor the effec- 
tive friction of a vibrated granular layer. If we choose the 
grain diameter d as the characteristic size and yd/g as 
the characteristic time, there remains three control pa- 
rameters, the amplitude a, the angular frequency u> and 
the cross section surface S of the droplet and one local 
parameter evolving in time, h/d. 

From the local slope of W(t) 3 , we have extracted a 
"local" effective friction coefficient and observed a sys- 
tematic decrease of /Lt(t), as H(t) decreases. Note that 
in the whole range of parameters, the coefficient /i re- 
mains within 30% of the dynamical friction coefficient 
fi d . Although a clear dependence on h is observed, the 
dynamics can be thought of as "adiabatic" , the shape 
remaining close to a parabola and the local slope of 
W(t) 3 giving the friction in the central part of the drop. 
Fig. ([3])b shows the effective friction coefficient extracted 
from spreading experiments at constant external driving 
parameters, starting from different initial shapes and dif- 
ferent surfaces S. A tiny systematic dependence on W 
may be evidenced, indicating limit of the adiabatic as- 
sumption. One could iterate the analytical procedure in 
order to include explicitly the local height dependence 
in the spreading law but this goes beyond the scope of 
the present paper. Hundreds of experimental measure- 
ments were performed that may all be described by the 
phenomenological formula: 



M 



Moo + (Mo - Moo) exp (-H/ {Md)) 



(9) 



treating the friction for very deep layers Moo as a f ree 
parameter, but with the same values for (1q ~ 0.37 and 
J\f ~ 7.8. Figure |3~l shows the fit coefficient 

Moo a s a 

function of the dimensionless parameter that collapses 
the best our data, I — -^=. The physical implications 
of this result is discussed later on. Using a parabolic ap- 
proximation of ^^(J), we extract for all the experiments 
the thickness dependence of the effective friction (fig.|3ji). 
We are finally able to propose an empirical friction law: 



M = Mo + (MP (I) - Mo) ({h/d) 



(10) 



that provides fair representation of our spreading exper- 
iments. 

Interpretation - We have designed an experiment that 
is, to our knowledge, the first explicit measurement of a 
granular rheology under controlled vibration. Similarly 
to sheared dense granular flows, the rheology depends on 
an inertia number I built as the ratio of a shear velocity 
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FIG. 3: (a) Function f3(T / [i) relating the sliding velocity to friction (see Eqs. ph,® and Q). It allows to measure an effective 
basal friction fi from the spreading dynamics, (b) Dependence of p on the drop height H for / = 24 Hz, V = 0.93 and different 
values of L x and S. The solid line is the best fit using eq. Q. (c) Effective friction fioo in the limit of a large height H as 
a function the rescaled vibration velocity I = acu/y/gd for / = 18 Hz (A), / = 22 Hz (V, / = 24 Hz (o), / = 26 Hz (o) and 
/ = 28 Hz (□). The solid line is the best fit by the parabola /ip(-T) = 0.67 - 0.22 (I - 1.5) 2 . (d) Height dependence of the 
effective friction extracted from the parabolic fit of all data: ((H/d) — (/ip — fj,) / (/ip — /io). Each square is an average over 
50 experiments. The solid line is the best fit by the form ((x) = 1 — exp (—x/Af), with M ~ 7.8. 



at the scale of the grain (here aui) and the typical impact 
velocity. However, this last parameter was found to scale 
on \J Pj p in dense granular flows [1]. Here, the presence 
of maximum in /ioo(-0 allows to dismiss this possibility. 
We find the impact velocity to scale on \fgd (a free fall 
from a grain size) instead of \fgh- Provided this surpris- 
ing difference, we observe a regime at low inertial number 
where the effective friction increases with /, as sheared 
dense flows [3]. Moreover, at larger /, a cross-over to a 
weakening regime is evidenced. Such an effect could be 
interpreted as a transition from a dense regime ('liquid 
phase) to a kinetic regime ('gaseous phase). Indeed, the 
generic prediction of kinetic theory is a decrease of the 
friction with the agitation (and thus /) [5]. The rheolog- 
ical form of Eq. (10), featured on Fig. [3J:, has recently 
been obtained from a mean field theory in |15j . In this 
model, the liquid regime is dominated by the trapping 
of grains and the gaseous regime by binary collisions, so 
that the collision rate increases at low I and decreases 
at large I. Another puzzling point is the finite size de- 
pendence when the layer gets thin. This nonlocal effect, 
involving about Af = 8 grains is reminiscent of the rhe- 
ology of thin sheared layers but here, goes qualitatively 
in the other direction: thinner layers bear less basal fric- 
tion. In further studies, it will be crucial to use this 'iner- 
tial tribometer' to investigate in details vibrated granular 
matter changing the nature and the shape of the grains 
and also getting closer to the jamming transition. 
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